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Abstract 

O I Geometries and energies for H3'(H2)^ clusters (n = 0, 11) have been cal- 



culated using standard "ab initio" methods. Up to clusters with n = 6, four 



I different Pople basis sets (DZ, TZ, TZP) have been used in the calculations. 

_ For larger cluster sizes, the calculations have been carried out with one basis 

cn 

CS| \ set (DZ) using the HF/CISD method. We discuss here only the nice coun- 



terplay of polarisation effects between the central Hg" ion and the adsorbed 



\ H2 molecules, which naturally governs both the geometric structure and the 

' energy of the clusters. 
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I. INTRODUCTION 



H3" molecules have been first experimentally verified in 1912 |T[|, Dawson and Tickner have 
identified Hj^(H2) clusters by mass spectometry in 1962 [Q] . In a recent work of Kirchner and 
Bowers on this topic 0] the dynamics of metastable H3 (Hg) fragmentations are studied. 
Another experimental work of Hiraoka and Mori [^] deals with the stability of H3'(H2)^ 
clusters. 

Many theoretical studies of small H3"(H2)„ clusters ( n < 3 ) have been published in 
the last years. Yamaguchi, Gaw, Remington and Schaefer |^ have made an intensive study 
of H3'(H2) with the Pople Double Zeta plus Polarization Basis (DZP), H^(H2)2 , H^(H2)g 
, H3" (112)4 ^^'-^ (112)5 have been calculated with the Triple Zeta plus Polarization Basis 
(TZP) by Farizon, Farizon-Mazuy, Castro Faria and Chermette 0, |§] and To our 
knowledge, this ab initio calculation is the first calculation for H3'(H2)^ clusters for n larger 
than 5. We will study here clusters up to n = 11 and compare the lighter ones the results 
of Farizon et al. The results will be discussed in terms of geometrical izomerization. 



II. NUMERICAL PROCEDURE 

The geometric optimization of the hydrogen clusters was implemented here by the 
BERNY-algorithm |TD[ . In the optimization no restrictions in bond-lengths and bond-angles 



have been adopted. The geometric construction of the clusters is done sequentially, starting 
with an optimized H3" structure with a randomly placed H2 molecule added at a distance of 
about 3.5 A from the . In a first run the position of the new H2 has been optimized, in 
a subsequent calculation all geometric variables have been optimized. The optimizations of 
the larger clusters are performed analogously with each new II2 molecule placed randomly 
at a distance of about 3.5 A from the last added H2 molecule. 

For the electronic configuration we selected the CISD/DZ and the CISD/TZP model for 
two sets of calculations. With regard to the electron-correlation effects in H2 and II3" , CI 
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calculations should be used [|TT] . 

For our intention to make predictions for larger clusters, we compare a model on a relatively 
low theoretical level but suited to be used for larger clusters, with a more sophisticated 
ansatz, both applied to small clusters. After comparing these results it should be possible 
to extrapolate properties of larger clusters using the simpler ansatz. 

For our extensive numerical calculations we adopted the GAUSSIAN 90 and GAUSSAIN 
92 code [jl2l on a S400/4 supercomputer |jl3|, which has a peak performance of 5 GFlops. 
Calculations of such relatively large systems with such large basis sets were only possible by 
use of an extra 1 GB ram-memory besides the main memory of up to 1 GB. The calculations 
took all in all about 200 cpu hours. 

The procedure described above has been tested by application to H3 and H3 (Hg) clusters, 
which have been investigated thoroughly before. 

The results of the calculations of the H3 (H2) cluster are in good aggreement with the results 
of Farizon et al. M] and with other theoretical work. 



III. RESULTS AND DISCUSSION 

The theoretical calculations, which include correlation effects of the electrons, reduce 
the energies of the systems. Regarding the electrostatic properties of the H2 molecules, it is 
nessesary to use polarization functions in these methods in geometric calculations. 

Fig. 1 and Fig. 2 display the electron density in H3"(H2) calculated with the TZP basis 
from different views. Apparently the H2 ist most tightly bound when closest to one of the 
hydrogens of pushing the negative charges to the other two hydrogens in the Hjj" . 
The more subtle effect of the orientation of the H2 is predicted by the TZP calculations 
to be perpendicular to the -plane. Thus the electrons of H2 are most distant from the 
electrons within the H3" . 

With the less sophisticated basis set DZ (two instead of six basis functions per electron) 
geometric optimization yields a quite different result. The polarization is less pronounced 
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and incomplete, yielding a completely planar structure. The electron charges in the Hjj" are 
almost equidistributed over the three protons. Interestingly the electron density between 
the H2 and the H;^ comes out too low, resulting in a reduction of the binding energy of as 
much as 0.05 a.u. . 

Understanding what happens in this smallest cluster makes it easy to understand the larger 
structures, too. 

The H3'(H2)2 and H3'(H2)3 clusters are built by the same structure principle, while the 
charge distribution is more amusing (the corresponding geometry structures are shown in 
figure 4 and figure 5). A nice example of charge frustration occurs in the , which means 
that the charges are being repelled from all H2 symmetrically and assemble in the middle of 
the H3" triangle. 

A first geometric and energetic shell is filled with these three molecules as is also indicated 
by the differences in the energies, which are give in figure ^ . 

For a study of the relative stability of the clusters we define a value which reflects that 
the H^(H2)„ clusters are a composition of H2 and molecules 

as a measure of the relative stability of the clusters. The relative stability reduces while the 
cluster size increases. 

Magic numbers using these defined /-values may be inferred for = 3, 5 and 9 using 
DZ. The relative stability of H^(H2)5 has been predicted by Hirao and Yamabe |14 
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geometric consideration. In the experimental work of Hiraoka ^ the II;^(Il2)g seems to be 
more stable than its neighbours. Hiraoka related this to a total planar geometry of H^(H2)g 
. We could not confirm this and found no planar structure. 

The magic numbers of n = 3,5,9, see figure 7, demonstrate the packing of the H2 , where 
three of them form the innermost shell, while the next shell consists of four II2 . 

In the calculations the model of the positive center of the cluster with neutal H2 
molecules has been verified. 



A property of the geometry of the clusters are the interatomic distances of the H2 molecules. 
This effect depends directly on the distance from the positive center of the cluster 
perturbated by the other H2 as shown in figure 6. The interatomic distances of the H2 
molecules decrease with increasing distance of those molecules from the H3 center. This 
is another nice polarization effect which can be understood in the following manner. The 
positive central charge attracts the electrons of the H2 reducing the screening of the protons 
of the H2 which in turn enlarges the repulsion between them. 

We have calculated H3"(H2)„ clusters up to n = 11, the largest H3"(H2)„ cluster calcu- 
lated so far. We observe two glooming obstacles: The optimization reveals more and more 
geometric isomers (second minima in the energy) with more and more shallow barriers. 
In addition the complexity of the basis set necessary for a reliable calculation should increase 
exponentially which is prevented by the computational resources. 

Finally the kinetic energies of the atoms contribute to the total energy of the clusters sub- 
stantially and are not taken into account in the codes used here. 
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FIGURES 

FIG. 1. Contour plot of the electron density of H3'(H2) with TZP basis taken in the plane of 

FIG. 2. Contour plot of the electron density of H3"(H2) with TZP basis taken perpendicular to 
the plane of H3" . 

FIG. 3. Contour plot of the electron density of H3"(H2) with DZ basis taken perpendicular to 
the plane of H3" . 

FIG. 4. Three dimensional view of the 113^(113)3 structure. 

FIG. 5. Three dimensional view of the H3" (112)3 structure. 

FIG. 6. H3"(H2)„ Geometric properties: interatomic distances of H2 as a function of the cluster 
size. 

FIG. 7. H3"(H2)„ Stabilities: 1- value of H3"(H2)„ clusters as a function of the cluster size. 
FIG. 8. Three dimensional view of the H3" (112)5 structure. 
FIG. 9. Three dimensional view of the H3" (112)7 structure. 
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TABLES 





TABLE I. Enorsios 


fa.u.) of n^flL, cliistor calculations 
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